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Abstract. We present a systematic survey of the capabilities of type la supernova explosion models starting from a number of 
flame seeds distributed around the center of the white dwarf star. To this end we greatly improved the resolution of the numerical 
simulations in the initial stages. This novel numerical approach facilitates a detailed study of multi-spot ignition scenarios with 
up to hundreds of ignition sparks. Two-dimensional simulations are shown to be inappropriate to study the effects of initial 
flame configurations. Based on a set of three-dimensional models, we conclude that multi-spot ignition scenarios may improve 
type la supernova models towards better agreement with observations. The achievable effect reaches a maximum at a limited 
number of flame ignition kernels as shown by the numerical models and corroborated by a simple dimensional analysis. 
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1. Introduction 

Over the past years, a consensus has emerged about the general 
astrophysical scenario of the majority of type la supernovae 
(SNe la). These events are associated with thermonuclear ex- 
plosions of white dwarf (WD) stars close to the Chandrasekhar 
mass (for a recent review see iHillebrandt & Niemeverl BoOO). 
A thermonuclear flame ignited near the WD's center is be- 
lieved to propagate outward incinerating most parts of it. 
The released nuclear energy may suffice to explode the star 
as could be s hown in several num e rical s i mulations (e.g. 
Nomoto et all 1 19841 iR einecke et al., l2002ct iGamezo et all 
2003URonke & HillebrandtH2005al) . 

Yet little is known about the way the thermonuclear flame 
ignites. The evolution towards flame ignition is a complex 
physical process. As a single WD is an inert object, dynamics 
must be introduced into the progenitor system by assuming it 
to be a binary. The favored scenario suggests a non-degenerate 
binary companion from which the WD accretes matter. Due to 
this mass accumulation it approaches the Chandrasekhar limit, 
steadily increasing its central density so that eventually carbon 
burning ignites. In the following several hundreds of years a 
simmering phase of convective burning sets the conditions un- 
der which finally a thermonuclear runaway occurs leading to 
the the formation of a flame. Unfortunately, the flame igni- 
tion is difficult to address both analytically and numerically. 
The few studies that tackled the pre-ignition evolution pro- 
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posed different scenarios for flame formation. While one study 
claims the flame ignitio n to take place in only one single point 
near the center fa6flich&Steirl E()02). others favor an igni- 
tion in multiple sparks distributed around the center or on only 
one side of it, depending on the large - scale convective flow 
pattern JCiarcia-Senz & Wooslevi 1 19951 IWooslev et all [2004; 
IWunsch & Wooslevll2004lllapichinoH2005h . 

Not being well constrained by theory, shape and location 
of the first flame(s) are usually treated as free initial param- 
eters in multi-dimensional simulations. The hope was orig- 
inally that those had only little impact on the models sup- 
porting the observational finding of a remarkable uniformity 
of SN la characteristics. However, multi-dimensional simula- 
tions starting with different initial flame configurations dis- 
agreed with this conjecture. The way of flame ignition turned 
out to be one of the influential parameters of the models. 
Different choices gave rise to controversial results. Single- 
spot central ignitions leadin g to exp l osions o f the WD were 
studied bv iNiem ever & Hillebrandl i 19951). IReinecke et al J 
Jl999ah. IHillebrandt et all J2000h. IReinecke et all J2002al) . 
IReinecke et al. J2002d) . and lGamezo et all J2003I) . In contrast, 
off-center single-s pot ignitions release insufficient energy t o 
explode the star jNiemever et all Il996l ICalder et all 12004 . 
Multi- sp ot ignition scen a rios wer e appli ed by iNjeineyereLalJ 
J 19961) IReinecke etaD J2002bl) . and iRopke & Hillebrandtl 
( 2005a) and seem to have the potential to increase the explosion 
strength. Apparently, a slight misalignment of the initial flame 
with the WD's center leads to a drastic change in the outcome 
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of th e simulations w hen starting with a single perfect sphere 
llCalder et all 12004). More structured multi-spot flame config- 
urations are less sensitive to it (Ropke & Hillebrandt, 2005a). 
But even setting aside the complications arising from asymme- 
tries in off-center ignitions, the results of the models have been 
shown t o depend on the numbe r of flame seeds jReinecke et all 
l2002bt ITravaglio etafl |2004 . The objective of the present 
study is to explore this effect in detail in a systematic approach. 

The question of what can be expected from multi-spot 
ignition scenarios with an increasing number and differ- 
ent distributi ons of initial flames wa s addressed only re- 
cently by IGarcfa-Senz & Bravol d2005l) . We report on a sim- 
ilar study, which is, however, based on a completely differ- 
ent approach to modeling t hermonuclear supernovae. While 
IGarcfa-Senz & Bravol J2005) apply a Lagrangian description of 
the hydrodynamics based on the smooth particle hydrodynam- 
ics (SPH) technique, we utilize an Eulerian grid-based finite 
volume method. The key distinction is that our approach facil- 
itates a self-consistent description of turbulent thermonuclear 
flame propagation, as will be discussed in Sect. [3] In contrast, 
the SPH model has to rely on a parameterization of the effective 
turbulent flame propagation velocity, since it cannot provide a 
valid description of turbulence effects. The predictive power 
of the SPH approach is thus limited. Nevertheless, the la rge- 
scale structures observed bv lGarcfa-Senz & Bravol |2005) ap- 
pear to be similar to what we find in our simulations, as can 
be seen from a comparison w ith a full-star model presented by 
Rop ke & Hillebrandn (12005 ah . This is not too surprising since 
both are driven by buoyancy and should be equally well repro- 
duced in both approaches. 

As will be discussed below, the simulation parameters 
and the details of the results of the present study dif - 
fer significantly from those of Gar cfa-Senz & Bravol J2005). 
Mor eover, the focuse s of the two surveys are distinct. While 
IGarcfa-Senz & Bravol d2005l) presented the first attempt to 
modeling SNe la in three dimensions (albeit in a parametrized 
way) from flame formation on and to assess the resulting explo- 
sions, we take a more pragmatic point of view. Our objective 
is to answer the question whether multi-spot ignition scenarios 
are in principle capable of curing some of the shortcomings of 
current deflagration SN la models. We therefore apply a num- 
ber of different initial flame configurations without attempting 
to model their pre-ignition evolution. The question of how re- 
alistic these are is set aside in the present study. 

Although being within the range of observational expecta- 
tions, the explosion energies and the masses of burning prod- 
ucts in all multi-dimensional deflagration models presented so 
far seem to be on the weak side and multi-spot ignition scenar- 
ios may offer a way to improve the results. Besides such global 
quantities, we are especially interested in the distribution of the 
species in velocity space. Previous work showed that poorly 
resolved models are inconsistent with observations in predict- 
ing low-velocity ox ygen and carbon lines in late time spectra 
(Koz maetalll2005l) . This is due to downdrafts of unburnt ma- 
terial in between burning buoyancy-driven bubbles of ashes. 
These downdrafts transport significant amounts of carbon and 
oxygen towards the center of the WD. It may well be, however, 
that the downdrafts carry sparks seeding additional burning, 



which are not resolved in current simulations. Moreover, the 
simu lation from which the synthetic spectrum of Ko zma et all 
(2005) was derived started out with a very artificial initial cen- 
tral flame configuration. The model was calculated on a uni- 
form computational grid of [256] 3 cells co-expanding with the 
WD. This led to an only marginally converged result and the 
explosion energy as well as the production of iron group and 
intermediate mass elements was rather low. Obviously, some of 
the disagreements with the observations are caused by the sim- 
plicity of the model - in particular the initial flame shape - and 
one may wonder whether the problems are mitigated with more 
structured flame configurations such as arising from multi-spot 
ignition. 

An estimate of the capabilities of the multi-spot ignition 
scenario is given in Sect.|2] In Sect.|3]we will briefly outline the 
techniques underlying our numerical simulations and discuss a 
new implementation that enables us to study multi-spot igni- 
tion scenarios in detail. This implementation is compared with 
previous simulations in Sect. |4] Results from two-dimensional 
models will be discussed in Sect. [5] Although it will be shown 
that these are inappropriate to study multi-spot ignitions, they 
point to important effects. A systematic survey based on three- 
dimensional simulations is presented in Sect. |6] In Sect.0we 
draw conclusions for modeling thermonuclear supernova ex- 
plosions. 

2. Multi-spot ignition: estimating the gains 

IWooslev etai1 ( l2004 and lWunsch & WooslevN2004 conclude 
from analytical models of the ignition process that a multi-spot 
ignition is possible (but not guaranteed). Timescale arguments 
show that a number of hot spots is in principle capable of evolv- 
ing towards a thermonuclear runaway. The first flame will ig- 
nite at a radius ~ 150 km off-center of the WD and more can 
follow. Yet their number and spatial distribution cannot be con- 
clusively constrained by the analytical models. 

Thus, it is justified to regard these as free parameters in 
simulations of the explosion process. In this spirit, we address 
the question of the number of ignition spots in the present 
study. Concerning the spatial distribution we simplify matters 
by assuming a spherically symmetric probability density of ig- 
nitions to occur around the WD's center. This ignores possi- 
ble large-scale anisotropies of the ignition process due to low- 
order modes in the c onvective flow pattern prior to ignition 
dWooslev etall botM). The impact of such effects is subject to 
a forthcoming publication. 

Given the ambivalence of the ignition conditions, it is a le- 
gitimate question to ask which multi-spot configuration would 
lead to an optimal fuel consumption producing the most vigor- 
ous explosion and potentially burning most of the fuel near the 
center. A dimensional analysis such as presented in the follow- 
ing can shed some light on this issue. Sects. |3]and|6]will tackle 
the question from the side of numerical simulations. 

One has to note that there are two major effects on the out- 
come of multi-spot ignitited supernova simulations. If a single 
flame seed is separated from the bulk of ignition sparks towards 
larger radii, it will experience a larger gravitational accelera- 
tion. The burning front resulting from this spark will evolve 
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faster than the other ignition points and dominate the flame evo- 
lution. This effect is determined by the nonlinear evolution of 
the burning fronts and hard to predict analytically. It will thus 
be discussed on the basis of numerical simulations in Sects. [5] 
and|5] Here, we only note that the distribution of ignition points 
should not be too sparse in order to avoid such effects. 

On the other hand, accomodating too many ignition kernels 
in a given volume will have the effect that the flames merge 
shortly after ignition due to self-propagation. In this case, the 
simulation will look similar to those ignited centrally in a sin- 
gle connected shape. The main advantage of multi-spot ignition 
scenarios, i.e. a large flame surface, is lost. Two effects counter- 
act the surface distruction by merging of the fronts. Firstly, the 
overall expansion of the star due to the nuclear energy release 
leads to an increasing separation of burning bubbles and, sec- 
ondly, buoyancy-induced flotation rises them in radial direction 
further increasing their separation. 

In a simplified picture, two igniting bubbles of radius r b at 
a separation / will lead to a maximum flame surface if their 
growth due to flame propagation is exactly compensated by 
these two effects. Assuming a homologous expansion of the 
star, the temporal change in distance of the bubbles, / ex induced 
by this effect is given by 



I - 2r b 



(i) 



where Td yn ~ {GpT 1 ^ 2 = 0.07 s is the dynamical timescale as- 
suming a density of 3 x 10 9 g ctrT 3 (G denotes Newton's con- 
stant). 

Due to buoyancy, two neighboring bubbles enclosing an an- 
gle <p with the star's center will rise radially, increasing then- 
separation by 



/buoy = 2rsin -. 



For the angle <p we assume 
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with an average initial distance (r) of the bubbles from the 
WD's center. Thus, 



buoy 



J ~ 2r b 

<r> 
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The asymptotic radial rise velocity of a bubble is (up to a factor 
of order unity) given by 



The distance decrement by growth of the bubbles due to 
burning is given by the propagation velocity of the flame front, 

'bum = 2Ub urn , (6) 

with i>bum ~ 10 7 cm s _1 shortly after ignition. Equating the sum 
of / ex and /buoy with /bum yields for the optimal initial bubble 
separation / opt 
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2UburnTdyn(>"} 
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The average density of bubbles n in a sphere of radius R is given 
by 

with N denoting the number of bubbles. In the optimal case, 
this should be equal to the inverse volume of a sphere with 
radius 0.5/ opt + r\,, i.e. the number of bubbles in this case is 
given by 



N pt = 



R ((r) + T dyn V 'At gr b ) 



TdynUbum<r> + 2r b <r) + 2T d y n r b V 'At g r b 



(9) 



Thus, assuming a central density of the WD of 3 x 10 9 ge m" 3 , 
an ignition radius of 100 to 150 km ( Woos lev et alJ. 12004). and 
setting the average initial distance from the center (r) to one 
half of this radius, the optimal number of ignition sparks with 
a radius of r b = 5.5 km (being the smallest resolvable bubble 
structure in the numerical models discussed below) should be 
somewhere between 23 and 76 per octant. This, of course, is no 
more than a crude estimate. Due to burning the bubbles will in- 
crease in size and this in turn affects the energy release and the 
contribution from buoyancy effects. Moreover, the establishing 
flows and the interaction of the flame with turbulence have been 
ignored. In a realistic numerical simulation, the flame seed con- 
figuration is much more complex. In most setups (see Sects. [5] 
and|6j we will apply a Gaussian distribution in radius. Since 
the outer bubbles are more separated in this case, the optimal 
number of bubbles is expected to be larger than the simple es- 
timate. 

With Eq. (|9} the question can be answered of whether it 
is possible to tune models to give more vigorous explosions 
by increasing the resolution and accommodating an ever larger 
number of ignition spots. Even in the limit of r b — » the opti- 
mal number 7Y opt of flame seeds increases to no more than 354 
to 1 195 per octant. Thus the effect of the number of initial flame 
sparks on the explosion strength is limited. 



(5) 



JPavies & Tavloilll95oT) . where At = (p a - p b )/(p u + Pb) de- 
notes the Atwood number across the flame front with the den- 
sity p u of unburnt material outside and the density pb of the 
burnt material inside the bubble. In its functional form, Eq. 
follows the expression obtained from balancing the buoyancy 
and the drag forces for a spherical bubble. Typical values are 
g ~ 10 9 cm s~ 2 and about 0.07 for the Atwood number at a fuel 
density of pb ~ 3 x 10 9 g cm~ 3 dTimmes & Wooslevifl992i) . 



3. Explosion: astrophysical scenario and 
numerical methods 

The explosion model describes the propagation of the flame 
after ignition near the center of the WD outwards. Our im- 
plementation is based on the deflagration scenario in which 
the flame velocity is sub-sonic and accelerated by the interac- 
tion with turbulence generated by generic instabilities (see e.g. 
iHillebrandt & Nie meYer[|2 00Cl). It follows the num erica l meth- 
ods developed bv iReinecke et al.l ([l£99b, 2002a) and Ropke 




Fig. 2. Snapshots of the evolution of a ci-model on a 256 2 cells hybrid grid. The logarithm of the density is color-coded and the 
contour marks the location of the flame front. 



(2005). For the details and extensive tests of the implementa- 
tion we refer to these works. The fundamental concept behind 
the approach is that of large eddy simulations (LES). Only the 
largest scales of the problem are directly resolved. The interac- 
tion of the flame with turbulent eddies on unresolved scales of 
the establishing turbulent ene rgy cascade is described by a sub - 
grid scale turbulence model JNiemever & HillebrandA IT995 ). 
The large separation of the flame width (of the order of a mil- 
limeter) from the resolved scales justifies the flame modeling 
as a sharp discontinuity separating the fuel from the ashes. To 
this end, the level-set techn ique is applied as described in detail 
bvlReineckeetakldlQQQhl). 

We apply t he sam e WD equation of state as used by 
Reinecke et al ] j2QP2a). In all simulations presented in the 
present paper we fix the initial composition to equal parts of 
carbon and oxygen throughout the star and construct the ini- 
tial WD assuming cold isothermal (T — 5 x 10 5 K) conditions 
and a central density of 2.9 x 10 9 gctrT 3 . The nuclear reac- 



tions are implemented in the simplified approach suggested by 
iReinecke et all J2002ah . Five species ( 12 C, 16 0, 24 Mg repre- 
senting the intermediate mass elements, 56 Ni as a representative 
of the iron group, and ff-particles) are followed 1 . At high fuel 
densities the material crossed by the flame is converted to nu- 
clear statistical equilibrium modeled as a mixture of "Ni" and 
ff-particles. Below fuel densities of 5.25xl0 7 g cirT 3 intermedi- 
ate mass elements are produced and once the fuel density drops 
below 10 7 g crrT 3 we stop burning. In the current implementa- 
tion electron captures are neglected. These could decrease the 
density in the ashes and therefore increase the density contrast 
over the flame. In the setups chosen for our simulations we ex- 
pect the effect on the dynamics of the explosion to be small, 
but significant changes may result at still higher central densi- 
ties than those applied here. 



1 We will set "Ni" and "Mg" in quotes henceforth to avoid confu- 
sion with the actual isotopes. 



F. K. Ropke et al,: Multi-spot ignition in type la supernova models 



6x 10' 



4x 10' 



2 x 10 s 




50 



100 150 200 

computational cells 



300 



Fig. 1. Evolution of the hybrid computational grid in a three- 
dimensional simulation. Each curve corresponds to a time pro- 
gression of At = 0.1 s, starting at the lowest curve with t — s. 
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Fig. 3. Flame resolution in different simulations: 256 2 cells hy- 
brid grid (solid), 256 2 cells uniform grid (dotted), 512 2 cells 
uniform grid (dashed), and 1024 2 cells uniform grid (dashed- 
dotted). 



Our novel approach that enables us to study the effects 
of initial flame configurations in detail is a computational 
grid different from previous implementations. While the first 
three-dimensional simul ations by Hillebra ndt et alJ J2000h and 
iReinecke et alJ J2002alblcl) were carried out on static Cartesian 
grid geometries with a fine-resolved uniform inner part and an 
exponentially growing grid spacing further out to capture parts 
of the expansion in the expl osion process, Ropke (2005) and 
Rop ke & Hillebrandtl J2005al) applied a moving uniform grid 
that tracked the expansion of the WD star. Here, we combine 
both approaches and use a moving grid that is composed of two 
nested sub-grids. The inner part {grid 1) again features a uni- 
form fine resolution and the outer grid cells (grid 2) grow ex- 
ponentially. Initially, both grids are moved individually. While 
grid 1 contains the flame and tracks its propagation, grid 2 fol- 
lows the WD expansion. Since in the course of the explosion 
the flame spreads out over the WD star, grid 1 must expand 
faster than grid 2. It is therefore possible to subsequently gather 
cells of grid 2 into grid 1 as soon as the grid spacings match. 
Eventually, grid 1 incorporates grid 2 and the full WD is cov- 
ered by a single uniform grid tracking its expansion. We will 
call this approach hybrid grid in the following. An example for 
the grid evolution in one of our simulations is given in Fig.^ 

Although conceptually simple, the hybrid grid has proven 
to be extremely useful and is much easier to handle than adap- 
tive mesh refinement strategies. Since the expansion and flame 
spread in the explosion is on average spherical for most igni- 
tions scenarios, it offers the possibility to maximally resolve 
the flame region with a given fixed number of computational 
cells. This leads to an improved modeling of the flame propa- 
gation and, in particular, provides the possibility to drastically 
improve the resolution in the central parts at the onset of the 
explosion. In this way it becomes feasible to resolve highly 
structured initial flame configurations. For studying multi-spot 
ignition scenarios, the number of initial flame kernels could be 



substantially increased. Whereas the previous static grid setups 
of one octant of the WD star could res olve up to ~30 bubble s 
in a simulation with 768 3 grid cells (cf. Travag lio et allEoO^l) . 
a hybrid-grid setup with 256 3 cells can accommodate several 
hundreds of ignition spots. 

All simulations performed in this study are carried out on 
only one spatial octant of the WD assu ming mirror symmetry 
to the other octants. As demonstrated by Ropke & Hillebrandt 
(l2005al) . this symmetry constraint does not suppress a develop- 
ment of potential low- wavenumber modes in the nonlinear flow 
pattern. Although these may in principle occur in convective 
phenomena, the short explosion time-scale in SNe la prevents 
them. Therefore restricting the simulations to only one octant 
does not miss physical effects and provides a valid approach to 
save computer time in our study. 

4. Testing the implementation 

Apart from a better resolution of the ignition, our approach 
concentrates resolution in the flame region also in the subse- 
quent stages of the explosion. It is therefore worthwhile to com- 
pare the results with previous simulations. In particular, we will 
point out some cha nges with resp ect to the expanding uniform 
grid simulations o f Ropke (2005). For this comparison, we ap- 
ply the standard test case of a single-octant simulation in which 
the flame is centrally ignited with a sinusoidal (2d) or toroidal 
(3d) perturbati on impo sed on it. According to the notation of 
IReinecke etaD ( l2002al) . we will refer to this flame setup as c3. 

Figure|5]shows snapshots from the evolution of a c3-model 
in two dimensions imposing cylindrical symmetry. Here, a 
256 2 cells hybrid grid was applied. The initial c3 flame is de- 
picted in the upper left panel of Fig.|2] Comparing these snap- 
shots to those of the previ ous uniform-grid implementations 
(cf. Fig. 4 of Ropke, 120051) . we note that the global structure 
resembles the simulations there carried out on a much larger 
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number of computational cells. These, however, possess finer 
sub-structures. The reason is clear from a comparison of the 
resolutions of the flame fronts in the different implementations 
provided by Fig. [3] Due to the moving grids the resolution 
changes with time in all cases. As expected, the concentration 
of computational cells in the flame part leads to a much better 
resolution of the flame in the first stages for the hybrid im- 
plementation on 256 2 cells than for the uniformly expanding 
256 2 cells grid. It is similar to the uniformly expanding 1024 2 
cells grid implementation. Therefore, the large-scale structures, 
seeded by the initial c3 flame perturbations evolve in a balanced 
way, whereas in the older implementation with 256 2 grid cells 
the inner Ray leigh-Taylor finger was suppressed (cf. Fig. 4 of 
lRopkeL E005). Since the flame propagates faster than the WD 
expands, the resolution in later stages - when the two nested 
grids have evolved into a single uniform grid - converges to- 
wards that of the previous implementation. Naturally, this pre- 
vents the flame from developing as detailed fine structures as 
the uniform grid 1024 2 implementation. This effect is, however, 
compensated with regard to the global quantities by the sub- 
grid scale turbulence model. We conclude that our novel imple- 
mentation provides a reasonable compromise between compu- 
tational expenses from larger grids and resolution of the flame. 

The situation is very similar in three-dimensional simula- 
tions. A snapshot from the standard c3_3d model 1.0 s after 
ignition is shown in Fig. 0] As in the two-dimensional case, 
the three-dimensional implementation of the hybrid grid leads 
to a more balanced evolution of the flame morphology im- 
printed on the initial flame seed. The c3_3d model in pr evious 
uniform grid implementations (see Fig. 13 of Ropke, 2005) 
showed three strong flame features evolving along the axes due 
to a suppression of initial perturbations. This grid-imprinted 
symmetry is relaxed in the novel implementation in which the 
flame is structured into more and smaller features (see Fig.|4j. 
Interestingly, the energy production is not increased by these 
effects in the three-dimensional simulations (cf. the energy re- 
lease of t he c3 3d mo del on a uniform expanding grid pre- 
sented by iRopkel Eo05 ). corroborating the conclusion that the 
description of the burning (at least in the stages where it pro- 
ceeds to iron group elements) is numerically converged in our 
models. Although increased resolutions give rise to a modi- 
fied structure of the flame and a redistribution of the ashes, the 
global quantities are largely unaffected. 

5. Two-dimensional simulations 

With the hybrid grid implementation, four different flame con- 
figurations were tested in two-dimensional simulations again 
assuming a cylindrical symmetry. All were carried out on 256 3 
grid cells with an initial inner grid spacing of 1 .74 km. Apart 
from the ci_2ii-model, different numbers of initial bubbles 
were chosen in the following procedure. The bubbles with a 
radius of 5.5 km were distributed randomly in angle and ac- 
cording to a Gaussian probability distribution with a dispersion 
of <t — 75 km in radius. This Gaussian radial distribution, how- 
ever, was distorted by the constraint that no bubble was allowed 
to ignite at a radius larger than 2.5cr and by imposing a mini- 
mum distance of the bubble centers af m j n . The models b30_2d, 
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Fig. 4. Flame front in model c3_3d 1.0 s after ignition. 



bl00_2d, and b200_2d contained 30 bubbles with d m [ n = 0.8, 
100 bubbles with d m [ a = 0.2, and 200 bubbles with c/ m ; n =0.1, 
respectively. 

The initial flame configurations are shown in Fig. [5] This 
figure also contains snapshots at t = 2 s. From a compari- 
son with the c3_2d-model in Fig.|2]it is clear that in the two- 
dimensional case multi-spot ignition scenarios do not give rise 
to an increased overall fuel consumption. This is confirmed by 
the evolution of the total energies in the different models plot- 
ted in Fig. [6] The models b30_2d and bl00_2d produce only 
about 40% of the asymptotic kinetic energy of model c3_3d 
and the energy production of model b200_2d falls in between. 

The reason for this finding is a peculiarity of two- 
dimensional models. Apart from the fact that turbulence in the 
two-dimensional case follows a different scaling than three- 
dimensional turbulence, two-dimensional simulations tend to 
amplify large structures. This may be attributed to the miss- 
ing degree of freedom in the third spatial direction where axial 
symmetry is imposed. Thus the initial flames consist of tori 
rather than bubbles. The large-scale flows around these struc- 
tures are expected to significantly differ from the flow patterns 
in three-dimensional simulations. 

As a result, in cases with few ignition spots, the outermost 
flame seed dominates the evolution. This is illustrated by snap- 
shots of the early flame evolution in model b30_2d in Fig. 
Starting from the initial configuration shown in Fig.|5j the inner 
flames merge quickly, but the outermost ignition spot (marked 
with the arrow in Fig. |5J evolves into a structure that domi- 
nates the flame evolution further on (cf. the snapshots in Fig.0 
and the snapshot at t — 2 s in Fig.|3. A similar effect governs 
the flame evolution in model bl00_2d. Due to the larger num- 
ber of ignition spots in model b200_2d the flame propagation 
is more balanced here and finally becomes dominated by two 
large features (see Fig. qj. The reason why the c3_2d-model 
appears so well-behaved is the symmetric and balanced pertur- 
bation we impose on the initial flame. It imprints three features 
which evolve equally. 
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Fig. 5. Two-dimensional multi-spot simulations: model b30_2d (top), model bl00_2d (middle), and model b200_2d (bottom). 
The left column shows the initial flame configuration and the right column provides snapshots of the evolution after 2 s with the 
the logarithm of the density color-coded and the flame front represented by the contour. The arrow in the top left panel marks the 
bubble that later on dominates the flame evolution. 
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Our results are in r easonable agreement with t he findings 
■ 1 j 1 c I n I 

of Niemever et al. ( 1996) and the recent study by Livne et al. 
(2005), both reporting on dominating large-scale flame fea- 
tures. 



Thus our findings point to the effect mentioned in Sect. |2] 
Single spots may lead to large-scale features which have 
an unfavorable effect on the burning. Since flow patterns in 
two-dimensional simulations differ significantly from those 
in three-dimensional cases, the question arises of how pro- 
nounced this effect is there. 



6. Three-dimensional simulations 

6.1. Ignition setups 

With the three-dimensional implementation of our scheme, we 
performed seven simulations with different initial flame con- 
figurations. The three-dimensional analog of our standard test 
model c3_2d, termed c3_3d, served as comparison with two- 
dimensional runs and with previous three-dimensional setups 
on other grid geometries. Five ignition configurations with bub- 
bles of a radius of 5.5 km were chosen in a similar way as the 
setups for the two-dimensional simulations of Sect. [5] Again, 
a Gaussian radial distribution was applied, with a dispersion 
of cr = 75 km. According to the number of igniting bubbles, 
the models are denoted as bl5 3d, b30 3d, b55 3d, bl50 3d, 
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^ 1 1 dimension is small. Therefore the shear flows are weak and 
the effects of turbulent flame wrinkling increase only gradually. 
This corresponds to a slow increase of the released energies as 
can be seen from Fig. [8] 
^ .„ At t ~ 0.4s the energy generation rates increase drasti- 
cally and peak around 0.6 s. The flame fronts at this time are 
shown in the middle column of Figs. l9l and 1 1 01 We note here 
that the initially separated ignition spots have mer ged into a 
single c onn ected structure in agreement with iReinecke et alJ 

- d2002bh andlRonke & H illebrandll d2005ah . As pointed out by 
, 2d Ropke & I lillehrandt (2005al) (see also Fig. 5 of that publica- 

b30 2d tion), this is a natural consequence of the flow field that es- 

bl00_2d tablishes around buoyantly rising burning bubbles. The stream 

b200_2d lines are directed around the forming mushroom-like struc- 

ture and converge behind it, where an upwards pointing flow 
j 5 ? Q emerges. This flow drags underlying flame patches towards 

the bubble so that the originally disconnected flame structures 
eventually merge. This effect has been ignored in the dimen- 
Fig.6. Total energies in the two-dimensional simulations. sional estimate of Sect l2l 

, , After t ~ 1 .5 s the burning has ceased in our models and the 

: -7= v = ;^ ;^_-JT^ energies released in the explosion process have reached their fi- 

- nal values (cf. Fig. [8}. The flame fronts at t — 2.0 s are shown 

- in the right columns of Figs. |9]and[n)] All models were fol- 
lowed up to t — 10.0 s where homologous expansion is reached 
with reasonable accuracy (Ropke, 2005). In the last seconds the 

~_ distributions of the ashes change only slowly in the relaxation 
; process. 

C3 f From Figs.l9landlT0lit is obvious that the different ignition 

b30 ~ 3d configurations lead to significant variations in the flame evolu- 

b55 3d ~ tions. They are determined by the effects mentioned in Sect. 13 

bl50_3d - The first effect is most relevant in cases of small numbers of 

b250_3d Z ignition spots. A sparse ignition may easily lead to anisotropics 

b500_3d - m tne evolving flame structure. Due to the randomness of flame 
j kernel locations the emerging large-scale flame structure is al- 
1.5 2.0 ready imprinted in the ignition seed. This is most obvious in 

model b55_3d, where a preferential direction towards observer 
Fig. 8. Total energies in the three-dimensional simulations. (cf. Fig|9] lower row) is already present in the ignition configu- 
ration and retained in the later evolution - although somewhat 
moderated in the latest stages. Since there are only very few ig- 
and b250_3d. In yet another setup, b500_3d, we assumed an nition spots at large radii, i.e. subject to the large gravitational 
equipartition over the radius instead of a Gaussian distribution, accelerations, these tend to produce the largest structures. This 
The models are summarized in Table [T] and the distributions effect is, however, less pronounced than in the two-dimensional 
of the ignition kernels in the multi-spot scenarios are depicted simulations. Thus, three-dimensional multi-spot models give 
in the left columns of Figs. |9] and [10| Potentially, the maxi- rise to a more robust flame evolution. 

mum distance of flame ignition is a relevant parameter since As expected, with increasing number of ignition spots the 
the gravitational acceleration increases steeply in the inner part flames will merge quickly due to burning and consequently 
of the WD. The values of the maximal ignition radii are given flame surface is lost. To illustrate this effect, the temporal evo- 
in Table[T] Fig.|8]shows the energetic evolution of the models. lution of an approximate measure of the flame surfaces in the 

different models is plotted in Fig. ^2 Obviously, they vary con- 
siderably. As a consequence of the dense distribution of initial 
flame bubbles in model b500_3d, the effect is most pronounced 
The evolutions of the flame fronts in our three-dimensional there. From ignition to t — 0.05 s the flame surface area in- 
models with multi-spot ignitions are shown in Figs. [9] and [H)] creases rapidly, but it again decreases to a minimum between 
As expected, the basic features are the same in all mod els and t = 0.1 s and 0.2 s. The flame surface destruction effect is also 
consistent with the findings of Reinec ke et al.l 1 12002b]) . After present in a less pronounced way in the models b250_3d and 
ignition in multiple spots, the volumes of the bubbles slowly bl50_3d and even b55_3d. In these, the flame surfaces increase 
increase due to the burning of the flame. In this first stage, monotonically, but the slopes decrease between t = 0.05 s and 
buoyancy-induced flotation of the bubbles is slow since their t = 0.15 s. In contrast, the surface destruction is not notice- 




L.0 

r[s] 



6.2. Flame front evolution 



10 F. K. Ropke et al.: Multi-spot ignition in type la supernova models 



Table 1. Flame configurations and results of three-dimensional models. 
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Fig. 11. Flame surface area in different models as a function of 
time. The inset enlarges the evolution in the early stages. 



able in simulations starting from sparse distributions of igni- 
tion spots. The slopes of flame surface area evolution increase 
monotonically in models bl5_3d and b30_3d. Interestingly, the 
c3_3d model corresponds to an intermediate case. 

The comparison of the later flame area evolution is more 
complicated. After t ~ 1.0 s the flame surface areas of the mod- 
els with sparse ignition seeds catch up with the models ignited 
in denser initial flame configurations. One possible explana- 
tion may be that in the latter, as soon as the initial flame seeds 
have merged due to burning, a new structure emerges which is 
similar to a single perturbed sphere. Since this large lump of 
ashes is generated in a random way, seeds for larger and per- 
haps dominating flame features may not emerge in a balanced 
distribution. Thus, the flame evolution may again be hampered 
by large scale modes as in the case of starting with sparse flame 
seeds and thus the loss in flame surface area may never be re- 
covered. 

It should be noted that the radius inside of which the flame 
kernels are ignited plays a minor role. In model b55_3d it was 
reduced by about 10% with little impact on the results. 



6.3. Global quantities 

The evolutions of the total energies in the models are plot- 
ted in Fig. [S] Values of the total energies, the nuclear energy 
releases, and the masses of produced iron group and inter- 
mediate mass nuclei and unburnt material at the end of our 
simulations (f = 10.0 s) are included in Table [l] The masses 
of unburnt material in all models are too large to be consis- 
tent with observations. This is partially due to the fact that 
we stop burning artificially once the fuel density drops below 
10 7 gem 3 . Realistically, burning should continue to lower fuel 
densities converting more carbon/oxygen material into inter- 
mediate mass el ements thereby increasing the energy release 
of the explosion (R opke & HillebrandtLl2005bT) . 

The most vigorous explosion resulted from model bl50_3d 
which released 1.18 x 10 51 erg of nuclear energy producing 
0.667 M of iron group elements and 0. 167 M of intermediate 
mass elements. The weakest explosion was obtained with the 
c3_3d initial flame configuration, but the model with most ig- 
nition spots, b500_3d, exploded almost as weakly. These global 
characteristics are a natural consequence of the flame area evo- 
lutions in the different models described above. The fact that 
model bl5_3d finally develops the largest flame area is no con- 
tradiction. It catches up with the flame area of model bl50_3d 
at t ~ 1.3 s, which is late in the burning phase of the simu- 
lation. At this time the burning is incomplete due to the low 
fuel densities in the expanded WD. It terminates in interme- 
diate mass nuclei and the energy release is thus lower than in 
complete burning to nuclear statistical equilibrium (NSE) in 
earlier stages. 

The global characteristics of the explosion models con- 
firm the conjecture of Sect. [2] The explosion strength one can 
achieve with multiple ignition spots is indeed limited and there 
exists an optimal flame configuration with a finite number of 
initial flame kernels. Nevertheless, most multi-spot ignition 
scenarios show a significantly increased energy production as 
compared to the centrally ignited c3_3d simulation. It is re- 
markable that - in contrast to the two-dimensional models - 
most of the simulations reach total energies within a narrow 
range of about 0.7 x 10 50 erg around 6.3 X 10 50 erg (cf. Fig.OD). 
This may be interpreted as a sign for convergence and robust- 
ness of the multi-spot ignition model in its three-dimensional 
implementation. As far as the global quantities are concerned, 
even the most energetic explosion model bl50_3d is not a 
clearly distinguished optimum, but all models starting with 15 
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Fig. 9. Evolution of the flame fronts in models bl5_3d (top), b30_3d (middle), and b55_3d (bottom), at times t — Os (left 
column), t - 0.6 s (middle column), and t - 2.0 s (right column). The indicated length scale applies to the respective column. 



to 250 ignition spots per octant reach similar explosion ener- 
gies with some scatter due to the random choice of initial flame 
locations. 

6.4. Distribution of elements in the ejecta 

Besides the global quantities as the production of energy, 
other quantities can be compared with observational findings to 
judge the validity of SN la explosion models. These are mainly 
the abundances and distributions of various species in the ex- 



plosion ejecta. Of course, most important is the radioactive iso- 
tope 56 Ni, since its decay powers the optical event, but also 
other isotopes contribute significantly to the shapes of spectra 
and light curves. 

Constraints on abundances and distributions of elements 
can be obtained from comparing synthetic light curves 
and spectra derived from the models with observations. 
Alternatively, one can determine the composition stratification 
of the ejecta by fi t ting m odel spectra to a series of observed 
ones (Stehle et al., 2005) and compare it with the results of 
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Fig. 10. Evolution of the flame fronts in models bl50_3d (top), b250_3d (middle), and b500_3d (bottom), at times t — s (left 
column), t = 0.6 s (middle column), and t = 2.0 s (right column). The indicated length scale applies to the respective column. 



explosion models. Such detailed approaches are, however, be- 
yond the scope of the present paper. 

Without nuclear postprocessing (cf. Trava glio et all 120041) 
we can only evaluate the cumulative abundances and distribu- 
tions of iron group elements, intermediate mass elements, and 
unburnt material (carbon and oxygen). The produced masses of 
the for mer are given in TableQ] From comparison with the find- 
ings of Trava glio et alJ d2004l) it seems likely that the amounts 
of 56 Ni synthesized in the most strongly exploding models will 
be close to 0.4 M Q . This value falls within the range of the ex- 



pectations for normal SNe la IContardo et allEoO O 1 ). but rather 
on the low side. 

The most interesting feature is the distributi on of un burnt 
materi al in velocity space plotted in Fig. El iKozma et alJ 
(2005) showed that spectra of the nebular phase put sharp con- 
straints on the models in this respect. A synthetic nebular spec- 
trum derived from a weakly exploding c3_3d model (but with 
a coarser flame resolution than the model presented here due 
to differences in the computational grids) featured strong emis- 
sion lines originating from oxygen at low velocities. This effect 
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Fig. 12. Distribution of the followed species in the three-dimensional models: carbon/oxygen fuel (solid) iron group elements 
(dashed), intermediate mass elements (dotted). 



might be an artifact of the peculiar (and unrealistic) c3_3d- 
configuration, where the large buoyant structures are already 
imprinted on the initial flame leading to strong downdrafts 
of unburnt material. It may be suspected that multi-spot ig- 
nition scenarios leading to different flame evolutions alleviate 
the downdraft problem. This could be possible due to the more 
complex initial flame structure which may be regarded as per- 
turbed by a wider range of modes. Moreover, due to the large 
flame surface area that is achievable here, burning becomes 
more efficient and increased amounts of the sinking fuel in the 
downdrafts may be converted. 

Indeed, this conjecture is supported by our results (cf. 
Fig. ll2> . In particular, in model b250_3d the unburnt material 
at low velocities is significantly reduced and the composition is 
dominated here by iron group elements. With increasing num- 
bers of ignition spots a clear trend is noticeable from Fig. 1121 
Sparse distributions of initial flames still show an unburnt- 
material dominated composition at low velocities. With in- 
creasing numbers of ignition spots this gradually changes into 
an iron-group dominance. This effect is, however, not simply 
caused by the fact that with denser distributions of flame ker- 
nels more material is converted already in the ignition setup. 
The trend weakens with larger numbers of initial flames (cf. 
models b250_c3 and b500_c3 in Fig. ll2l> . The reason is most 
likely again the effect of flame kernels melting together shortly 
after ignition, forming a large lump as described above. The 
emerging random perturbations of this structure are likely to 
prevent a balanced flame evolution similar to what is found in 
case of sparse ignition spot distributions. Such flame propaga- 



tion scenarios naturally leave behind more unburnt material in 
central regions. 

In contrast to the global quantities discussed in Sect. 16.31 
the distribution of the species in the explosion ejecta clearly 
depends on the number of flame ignition kernels. Models with 
~150 ignition spots per octant produce distributions that are 
closest to observational expectations. 

7. Conclusions 

In the present study we have explored the possibilities of multi- 
spot ignition scenarios of Type la supernova explosion models 
in a systematic way. Owing to improvements of the computa- 
tional code, it was possible to achieve a high initial resolution 
of 1 .74 km on computational grids with 256 cells per dimen- 
sion. By setting the initial bubble radii to 5.5 km, it was possi- 
ble to accommodate up to 500 flame kernels within a radius of 
~ 180 km around the center of the WD star. Due to the random 
placement and overlap of the bubbles, larger structures emerge 
providing a variety of sizes of ignition seeds. 

Comparing the results of two-dimensional simulations with 
those of three-dimensional models, we find significant differ- 
ences in the flame front evolutions starting from multiple igni- 
tion spots. In the former setups, single outliers tend to dominate 
the evolution to a much higher degree than in three dimensions. 
We conclude from these findings that two-dimensional simula- 
tions are inappropriate to explore initial flame configurations. 
Two-dimensional studies can therefore only be a tool to study 
other effects and compare simulations with the same fixed and 
carefully chosen initial flame configuration, a possibility be- 
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ing the c3_2d-model. Thus, although the general picture of 
our two-dimensio nal simulations is similar to the findings of 
iLivne et al.l (Eo05). we disagree with their conclusions. 

A simple estimate showed that the gain in explosion 
strength arising from multi-spot ignition scenarios should be 
limited. This is confirmed by a systematic study based on 
three-dimensional simulations covering one octant of the star. 
Different ignition configurations with 15 up to 500 flame seeds 
were tested. We found the most vigorously exploding model 
for 150 flame kernels per octant with a radial Gaussian dis- 
tribution inside a sphere of 180 km around the center of the 
WD. This agrees reasonably with our dimensional estimate in 
Sect. |2j given the simplifications made there. Setups with in- 
creased resolutions - although allowing for smaller flame seeds 
- would not profit from much higher numbers of ignition spots. 

In terms of production of energy and iron group elements 
this simulation came close to the feiO-model 2 o flTravaglio e t al. 
J2004I) . which is a success of the hybrid grid implementation 
since it allows results on an only [256] 3 cells computational 
grid that previously required [768] 3 cells. This novel approach 
will provide a powerful tool to study the details of SN la defla- 
gration models at moderate computational expenses. The sim- 
ilarity with the fe30-model suggests that the deflagration sce- 
nario of SNe la in its current implementation (i.e. applying a 
multi-spot ignition with a Gaussian radial distribution within 
~ 180 km around the center of the WD and ignoring burning 
below fuel densities of 10 7 gcnT 3 ) will not allow asymptotic 
kinetic energies much higher than ~7 x 10 50 erg and the to- 
tal production of iron group elements will not greatly exceed 
0.7 M Q . Since these values already fall in the rang e of obser- 
vational expectations (cf. eg. Contard o et allEoOOl) . multi-spot 
ignition scenarios offer a way of modeling "normal" SNe la 
events. These global quantities seem to be rather robust against 
variations in the number and configurati on of flame kernels. 

Comparing our results with those of lGarcfa-Senz & Bravo 
( 120051) we note agreement in some general results. Multi-spot 
ignition scenarios may give rise to an explosion of the WD star. 
The evolution of the large-scale features is similar in both stud- 
ies. Our novel approach, however, facilitated an increased res- 
olution of the initial flame configuration by a factor of 10 com- 
pared with Ga rcfa-Senz & Bra vo (2005j). It was therefore pos- 
sible to accommodate much smaller flame kernels in greater 
numbers within the ignition region. This may be the main rea- 
son for disagreements in the conclusions (our more elaborate 
flame description certainly contri butes as well, but this is diffi- 
cult to disentangle). In contrast to lGarcfa-Senz & Bravol J2005) 
our results do not indicate a convergence with ever increasing 
numbers of flame seeds, but rather a maximum at a limited 
number. Moreover, our initial flame bubbles seem to merge 
more rapidly forming a connected flame structure. This gives 
rise to differences in the distribution of species in the ejecta. 
The dumpiness of ashes reported by iGarcfa-Senz & Bravol 
J2005h is less pronounced in our simulations. 

In cases with larger numbers of ignition spots, we find a 
much better exhaustion of the fuel in central regions. In contrast 

2 Not to be confused with the b30_2d and b30_2d of the present 
study. 



to the global quantities, this effect is sensitive to the number 
of ignition spots, indicating an incompleteness of the current 
deflagration model of SNe la. Potentiall y, an improved model- 
ing of late phases of the burning process ( Ropke & Hillebrandd 
2005b) may be necessary, which possibly also increases the ex- 
plosion energy. 

We conclude from our simulations that multi-spot ignition 
scenarios may allow SN la models to achieve better agreement 
with observations, although the effect is limited with regard to 
the number of igniting flames. It thus does not lead to an ar- 
bitrary means of tuning the models to higher energy release or 
nuclear conversion in the explosive burning. Multi-spot igni- 
tion scenarios, however, may be capable of reproducing "nor- 
mal" SN la explosions. 
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